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Abstract. This paper reports about the development of two provably 
correct approximate algorithms which calculate the Euclidean shortest 
path (ESP) within a given cube-curve with arbitrary accuracy, defined 
by £ > 0, and in time complexity K,{e) ■ 0{n), where K{e) is the length 
difference between the path used for initialization and the minimum- 
length path, divided by e. A run-time diagram also illustrates this linear- 
time behavior of the implemented ESP algorithm. 

1 Introduction 

Euclidean shortest path (ESP) problems are defined by a (2D, 3D, ...) Euclidean 
space which contains (closed) polyhedral obstacles; the task is to compute a path 
which connects two given points in the space such that it does not intersect the 
interior of any obstacle, and it is of minimum Euclidean length. 

Examples are the ESP inside of a simple polygon, on the surface of a con- 
vex poly tope, or inside of a simply-connected polyhedron, or problems such as 
touring polygons, parts cutting, safari or zookeeper, or the watchman route. All- 
together, this defines a class of immensely important computational problems of 
huge impact in economy, science or technology. 

For time complexities of algorithms in this area, we cite two examples. The 
general 3D ESP problem (e.g., path- planning in robotics) is NP-hard, see J. 
Canny and J. H. Reif [5J. 

For 2D ESP problems, there are linear-time, but very complicated algorithms 
(e.g., algorithms for ESP calculation in a simple polygon, based on B. Chazelle's 
[6] triangulation of whole polygons), or linear-time and easy-to-implement algo- 
rithms (e.g., for the relative convex hull in the 2D grid, see [9]). 

In this paper we consider ESPs in simple cube-curves, which are formed 
by successively face-adjacent grid cubes (of the uniform orthogonal 3D grid, see 
digital geometry [TT]). T. Biilow and R. Klette published between 2000 and 2002 
(see, e.g., [4j) a so-called rubberband algorithm (RBA) for the calculation of a 
Euchdean shortest path in a simple cube-curve. [4] stated two open problems: is 
this approximate RBA actually always converging (with numbers of iterations) 
to the correct ESP, and is its time complexity actually linear as all experiments 
indicated at that time. 

This paper reports about the development of two approximate RBAs, which 
always converge towards the ESP, and have K{e) ■ 0{n) time complexity. The 



paper provides a first summary of work done by F. Li and R. Klette in 2003 
to 2006; for details and generalizations of the RBA approach for solving various 
ESP problems in 2D or 3D space, see |12ll3j and forthcoming publications of 
the authors. 

2 The Original RBA 

Critical edges of a given cube-curve g are those grid edges which are incident 
with three cubes of the curve (see Figure [ij . Critical edges are the only possible 
locations for vertices of an ESP [TU]. A subset of those will define the step set of 
the RBA, which contains all those critical edges which contain exactly one ESP 
vertex each. 




Fig. 1. Critical edges ei, 62, 63, 64, 65, and eg. 

The Original RBA, as published in |4lllj . is as follows: it consists of two 
subprocesses, (i) an initialization process (e.g., from an endpoint of one critical 
edge to the closest endpoint of the subsequent critical edge; satisfying a "closed- 
path" constraint at the end), and (ii) an iterative process which contracts the 
path during each of its loops, using a break-off criterion 

where e > 0, and £„ is the total length of the path after the nth loop. 

During each loop, the algorithm tries to shorten the path locally by checking 
three options, called OPl, OP2, and OP3. OPl and OP2 find the step set of 
critical edges. OPS optimizes the position of a vertex on its critical edge. These 
options are defined as follows: 

OPl: delete vertex pi if the line segment Pi-ipi -I- 1 is in the tube g, which 
is the union of all the grid cubes in the given simple cube-curve g; 

OP2: calculate intersection points of the triangle pi-iPiPi+i with all critical 
edges ("between" pi-i and Pi+i) and replace the subsequence pi^i,pi,pij^i by 
the resulting convex arc, defined by these of intersection points; 



OPS: move pi on its critical edge e into the optimum position Pnew, with 
de{Pnew,Pi-i) + de{Pi+i , Pnew) = mf {4 (p, Pi- 1 ) + 4(p-i+i,p) : P G e}, where de 
denotes the Euchdean distance. 

We continue with vertices Pnew,Pi+i,Pi+2 of the path. At the end of each 
loop we compare the total length of the new path with that of the path at the 
end of the previous loop. 

See Figure [2] for OP2. Here, vertices on critical edges en, ei4 and eig are 
replaced by a convex arc with vertices on critical edges en, 613, eig, and eis, 
and (in general) it may be en, 614 and eig again within a subsequent loop - of 
course, for a reduced length of the calculated path at this stage. 




Fig. 2. Illustration for the (original) Option 2. 



The situation with the original RBA in 2002 [1] was as follows: Even for 
very small values of e, the measured time complexity indicated 0{n), where n is 
the number of cubes in g. However, there was no proof for the asymptotic time 
complexity of the original RBA. For a small number of test examples, calculated 
paths seemed (!) to converge against the ESP. However, no implemented algo- 
rithm for calculating the correct ESP was available, and (more general) no proof 
whether the path, provided by the original RBA, converges towards the ESP. 
Nevertheless, the algorithm is in use since 2002 (e.g., in DNA research). 

3 Non-Existence of an Exact Arithmetic Algorithm 

An arithmetic algorithm consists of a finite number of steps of arithmetic oper- 
ations, possibly also using input parameters from the field of rational numbers, 
using only the following basic operators: +, — , •, / or the fcth root, for k >2. 

OPS can be formalized by a system of three PDEs, involving parameters 
ti £ R for critical edges e^ of the step set. The result ensures that pi{ti) is the 
optimum point on e^ . Considering the situation illustrated in Figure [3] this is 
equivalent to the problem of finding the roots of p{x) = 84x^ — 228a;^ -I- 361a;^ -t- 
20a;^ + 210a;^ -I- 200^ -|- 25. In fact, this problem is not solvable by radicals over 




Fig. 3. Consider the calculation of ti and t2 such that the polyline 
Po{to)pi{ti)p2{t2)p3.{h) is fully contained in g. 

the field of rationals; see [12]. (The proof uses a theorem by C. Bajaj [2] and the 
factorization algorithm by E. R. Berlekamp [3].) 

This example allows two corollaries. Obviously and well-known is that there 
is no exact arithmetic algorithm for calculating the roots of polynomials (known 
since E. Galois; B.L. van der Waerdens famous example is p{x) — ~ x ~ 1). 
And secondly, there is also no exact arithmetic algorithm for calculating 3D 
ESPs. C. Bajaj [T] showed this based on a polynomial of order 20 for the general 
3D ESP problem. As a new result, here we have an oder 6 polynomial, and the 
restricted ESP problem for simple cube-curves! 

Note that this is not just a "rounding number problem" but a fundamental 
non-existence of exact algorithms, no matter what kind of time-complexity is 
allowed. 

There is a uniquely defined shortest path, which passes through subsequent 
line segments ei, 62, . . . , in 3D space in this order; see, for example, [7j. Ob- 
viously, vertices of a shortest path can be at real division points, and even at 
those which cannot be represented by radicals over the field of rationals. 

4 Approximate Algorithms 

An algorithm is an (1 -I- e)- approximation algorithm for a minimization problem 
P iff, for each input instance / of P, the algorithm delivers a solution that is at 
most (1 -f e) times the optimum solution 0. 

The general 3D ESP problem can be solved in O {ri^ \b + log(7i/e)]^ /e^^ time 
by an (1 -f e)-approximation algorithm; see C. H. Papadimitriou |15j . 

An algorithm is K-linear iff its time complexity is in K(e) • ©(n), and function 
K does not depend on the problem size n, for e > 0. We use K{e) = (£0 — 
where C is the true length of the ESP, and Cq the initial length. 

A cube-curve is first-class iff each critical edge contains one ESP vertex. The 
original RBA is correct and K-linear for first-class cube-curves [12]. 



Fig. 4. Weighted undirected graph for m = 3. 



We analyzed the following approximate graph-theoretical algorithm: Subdi- 
vide each critical edge by m uniformly-spaced vertices; connect each vertex with 
those vertices such that the resulting edge is contained in the tube g. This defines 
a weighted undirected graph (see Figure |4| . Calculate a shortest-length cycle, 
and use this as a (first-class !) input for the original RBA. 

The time-complexity of the graph-theoretic algorithm (in our specification) 
equals O {m^ri^ + K{e) ■ nj. It applies Dijkstras algorithm repeatedly; possibly 
its time-complexity can be reduced, but certainly not to be K-linear. 

However, this (slow) algorithm allowed for the first time to evaluate results 
obtained by the original RBA. 

Assume a simple cube-curve g and a triple of consecutive critical edges ei, 
62, and 63 such that is orthogonal to ej, for i,j = 1,2,3 and i ^ j. If ei and 
63 are also coplanar, then we say that ei, 62, and 63 form an end angle, and a 
middle angle otherwise. 

The following approximate numerical algorithm requires an input which is 
first-class and has at least one end angle; the cube-curve is split at end angles 
into one or several arcs. For each arc, one vertex on each critical edge can be 
calculated using the systems of PDEs briefly mentioned already above; variable 
ti determines the position of vertex pi on edge e^. This algorithm is provably 
correct and K-linear for the assumed inputs. 

An open problem in |TT] (page 406) was stated as follows: Is there a simple 
cube-curve such that none of the vertices of its ESP is a grid vertex? The answer 
is "yes" [T^], and any of those curves does not have any end angle; see Figure 5p" 



Thus, the provably correct approximate numerical algorithm cannot be used in 
general. 

This lead us back to the initial two questions about the original RBA: is it 
correct? (We can use either the approximate graph- theoretical or the numer- 



Here are two new open problems: What is the smallest (say, in number of cubes or 
in number of critical edges - both is equivalent) simple cube curve which does not 
have any end angle? What is the smallest (say, in number of cubes or in number of 
critical edges - both is equivalent) simple cube curve which does not have any of its 
MLP vertices at a grid point location? We assume that the second problem is more 
difficult to solve. 



Fig. 5. A simple cube-curve where the ESP does not have any grid-point vertex 
(and which has no end angle). 

ical algorithm for evaluation.) What is its time-complexity in general? Indeed, 
corrections were in place: 

OP2: if intersecting with the triangle pi-iPiPi+i and using the convex arc 
only, we may miss edges of the step set (see Figure [6] for such a situation) - more 
tests are needed, and this option was totally reformulated (for details, see [12] - 
the specifications require some technical preparations which cannot be given in 
this short paper). 

OPS: the vertex Pnew, found by optimization, may specify edges Pi-iPnew 
and PnewPi+i such that one or both of them are not fully contained in the tube 
of the curve; an additional test is needed (a simple correction). 




Fig. 6. A situation where the original Option 2 fails. 



1DD 2D0 300 400 5D0 BDO 

X(total rumber of cubes) 



Fig. 7. Edge-based RBA Implemented in Java, ran under Matlab 7.0.4, Pentium 
4, using e = 10"^^°. 



Thus, as an end to this story right now, those corrections define a provably 
correct (for any simple cube-curve) and k- linear edge-based RBA 12J. 

Instead of moving points along critical edges, we can also move points within 
critical faces (which contain one critical edge). Of course, the vertices will finally 
move onto or towards critical edges. This conceptually simpler (in its OP2) face- 
based RBA is also provably correct, but showing a slower convergence (within 
the limits of being K-linear) towards the EPS. 

See Figure [T] for some statistics about measured run time. Half of a simple 
cube-curves was generated randomly, and the second half then generated using 
three straight arcs for closing the curve. The number of cubes in generated curves 
was between 10 and 630. The break-off criterion was defined by e = 10" 

Figure |8] illustrates the meaning of the break-off criterion. The lengths £„, 
for loops n = 1,2,3 .. . define a Cauchy sequence which converges towards the 

' ' Length 




Iteration 



Fig. 8. Let £ be the maximum accuracy of the program, that means the smallest 
number for discriminating between £„ and Cn+i - Still, the difference to the true 
value C might he S > e. The algorithm allows to obtain arbitrary accuracy (with 
respect to £) when continuing iterations, but this would require to reduce e. 



true length C. An in-depth study of this sequence may reveal whether we can 
assume 5 < e in general, or not. 

5 Conclusions 

This paper reported about the process of solving one particular ESP problem. 
The developed methodology [i.e., define "critical" subsets, specify the step set 
such that each critical subset in this set contains exactly one (possibly redun- 
dant, such as colinear) vertex, apply OPS] can be applied to ESP problems as 
considered (e.g.) in [TT. A few RBA applications have been illustrated in |12ll3j . 
and further applications will be published soon by the authors. 
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